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Abstract 

Hamilton variational principle for a special type of statistical ensemble of 
deterministic dynamical systems is derived. This form of variational principle 
allows one to describe the statistical ensemble in terms of wave functions 
and provides a basis for a description of statistical ensemble of stochastic 
systems. It is shown that sometimes such a statistical description of stochastic 
particle motion appears to coincide with the quantum description in terms of 
Schrodinger equation. 

PACS numbers: 03.40. Gc; 47.10.+g; 03.65.Bz 



1 



I Introduction 



One derives the Hamilton variational principle, which appears to be an effective 
mathematical tool of a statistical description. It allows one to write dynamic equa- 
tions for a statistical ensemble in terms of wave functions and to show that the 
quantum mechanical description in terms of the Schrodinger equation is a special 
case of statistical description. 

Let a deterministic dynamic system iSdQ be described by the Lagrangian L(t, x, dx./ dt) , 
where x = {x a } a = 1,2, ...n, are generalized coordinates and v = dx/dt = 
{dx a /dt} a = 1,2, ...n is generalized velocity. For brevity we shall speak about 
a particle with a position x and velocity v. A statistical ensemble £ gen [S&] (i.e. 
a set of many independent identical systems Sd) is a deterministic dynamic sys- 
tem, whose state is described by the distribution function F (t, x, p), where p = 
{dL/d(dx a /dt)} , a = 1,2, ...n is the generalized momentum of The state 
F (t, x, p) of the general ensemble S gen [«Sd] evolves according to the Liouville equa- 
tion 

dF_ _d_ f^_ F \ ( dH 

dt dx a \dp a J dp a \dx° 

where H = H (t, x, p) = vp — L is the Hamilton function for the dynamic system Sd- 
Summation is produced over repeating Greek indices (1 — n). Dynamic systems Sd, 
constituting the statistical ensemble S gen [Sd] , are called elements of this ensemble. 



£ S en[S d ): — + —[—F)- — [—F)=0 (1.1) 



Dynamic equation ( |1 . 1| ) for the statistical ensemble 8 gen [Sd] is obtained formally 
from the dynamic equations for elements Sd of the statistical ensemble S gen [Sd] ■ 

We are interested in a special type £ v [Sd] of the general ensemble S gen [Sd] which 
will be referred to as a pure statistical ensemble, because its state may be described 
by the wave function^. By definition a pure statistical ensemble £ gen [Sd] of dynam- 
ical systems Sd is such an ensemble, whose state F p (t, x, p) may be represented in 
the form 

S p [S d ] : F p (t, x, p) = p (t, x) 5 (p - P (t, x)) (1.2) 

where p (t, x) and P (i, x) = {P a (t, x)}, a — 1, 2, ...n are functions of only time t and 
generalized coordinates x. In other words, the pure ensemble £ p [Sd] is a dynamic 
system, considered in the configuration space of coordinates x. 
Dynamic equations for the pure ensemble S p [Sd] have the form 

l^t'ff^H (L3) 

d d ( dH \ dH 

W^ + ft? { pP °8R. (t ^ ¥) + "S^ (i ' X ' P) = ' /3= 1,2,.. .n (1.4) 



lr The term "dynamic system" is used as a collective concept with respect to terms "nondeter- 
ministic (stochastic) dynamic system" and "deterministic dynamic system". Such an unusual temi- 
nology is used, because we need a collective concept for concepts "dynamic system' and "stochastic 
system" . 

2 In quantum mechanics a pure state of a quantum system is called a such one which can be 
described by means of wave function (not by the dense matrix) jl| . 
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They are obtained after substituting (jl.2|) into (|1.1| ) and integration with respect to 
p with the multipliers l,pp, (3 = 1,2, ...n respectively. Interpreting p as a particle 
density and v = OH/ dp as a generalized velocity, the equation (|1.3|) is regarded as a 
continuity equation. Then ( |1.4|) may be interpreted as generalized Euler equations 
for some ideal fluid without pressure. 

Use of pure statistical ensembles is convenient in the case, when elements of the 
statistical ensemble are nondeterministic (stochastic) dynamic systems S st , for which 
there are no dynamic equations. Then the following procedure is used. The general 
statistical ensemble £ gen [S st ] is considered to be consisting of elements £ p [S st ], which 
consist in turn of elements S st , i.e. 

£gen [<Sst] — £gcn [£p [<S s t]] (1-5) 

There are no dynamic equations for dynamic system S st , but there are dynamic 
equations for a statistical ensemble £ p [S st \. Derivation of dynamic equations for 
£p [Sst] is rather complicated informal procedure, which needs a special considera- 
tion for any kind of stochastic system S st . But if dynamic equations for dynamic 
systems £ p [S st ] have been derived, dynamic equations for the statistical ensemble 
£gen [£ P [<Sst]] of deterministic dynamic systems £ p [S st ] are obtained as a result of 
a formal procedure. Thus, complicated informal procedure of dynamic equations 
derivation for S gen [S st ] is divided into two parts: (1) derivation of dynamic equa- 
tions for a more simple statistical ensemble £ p [S st ] by means of informal procedure, 

(2) derivation of dynamic equations for £ gen [£ p [S st ]] by means of a formal procedure. 
In this sense, a consideration of more simple pure statistical ensembles £ p sim- 
plifies investigation of stochastic systems. 

Let us add to the Euler equations the following equations 

dr 13 8H 

_ = _ (t , x ,P), , 8 = 1 , 2 ,...„ (1 . 6) 

describing motion of a particle in a given velocity field v =dH/dP. These equations 
can be rewritten in the form, known in hydrodynamics as Lin constraints 

^ + |^(t,x,P)cU /3 = 0, p = l,2,...n (1.7) 

where £ (t, x) = {£ Q (t, x)} , a = 1, 2, ...n are n independent integrals of n equations 

(03) • ' 

The system of 2n + 1 equations ( |1.3| ), (p..4|), ( |1.7| ) forms a complete system of 
dynamic equations, describing evolution of dynamic systems <Sd, constituting the 
pure ensemble £ p [S&\ . 

One can show that the system of 2n + 1 equations (|1.3|), (|L~4|), ( |1.7| ) reduces to 
the form of a system of n + 2 equations 

M<V + (f(£H£ Q ]+#(x,P) = 0, d ^ = ^ fc = 0,l,...n (1.8) 



3 



(dH \ 
p—(t,x,P)j =0 (1.9) 

-jf + gpfr x > p ) 9 °£f> = °> /5 = 1,2, ...n (1.10) 

where 93 is a new variable, and P is expressed now via n arbitrary functions g (|) = 
{g a (|)} , a = 1, 2, ...n of argument £. 

P,3 = 6 {dptp + (7° (<■)) fy£ a , /? = 1, 2, ...n (1.11) 

Practically it means that the system of 2n + 1 equations ( |1.9j ) , ( |1.4| ) , ( |1.10| ) is inte- 
grated partially. Here bo is an arbitrary constant, which can be incorporated in the 
variable (p and arbitrary functions g (£) . 

The system of n+2 equations (|1.8| ), (|1.9|) , ( |1.10| ) is remarkable in the relation that 



it can be described in terms of ip function (wave function). But it is more convenient 
one to transform these equations, using the Hamilton variational principle. 

We shall show that dynamic equations for the pure statistical ensemble £ p [Sd] 
of deterministic dynamic systems Sd are derived from the variational principle with 
the action 

£ p [<S d ] : A E \p,ip,Z\ = J p{-H(t,x,p) ~b [do^ + g a U)doU}d n+1 x, (1.12) 

Pp = bo[dp<p + g a {t)dpZ a ], /3=l,2,...,n,, d { = A (1.13) 

where p, ip, | are dependent variables, which are considered to be functions of x = 
{x°,x} = {t,x}. if(t,x, p) is the Hamiltonian function of Sd- &o is an arbitrary 
constant and g a {$,), a = 1,2, ...,n are arbitrary functions of the argument £. The 
dynamic variables ip, | are hydrodynamic (Clebsch) potentials. Clebsch || |] had 
introduced them for a description of incompressible fluid. The variables ip, £ are 
called potentials, because momentum p = P (t, x) is expressed via derivatives of 
potentials (p, £, as one can see from the relations ( |1.13| ). The Hamiltonian H (t, x, p) 
is a function which determines the form of the action (|1.12|) , and the variational 



principle based on the action ( |1.12j ) may be called the Hamilton variational principle. 

The dynamic system £ p [Sd] is an ideal fluidlike continuous medium without 
pressure. Variational principle for such a system in the Lagrangian coordinates £, 
(where £ are considered to be independent variables) is derived very simple. But 
a derivation of a variational principle for an ideal fluid in Euler coordinates x was 
being a problem for many years 0, because one tried to obtain a variational principle 
only for the Euler equations, i.e. for equations of the type (|L~3l) , ( |1.4|) . The Euler 
equations form a closed subsystem of the complete system of dynamic equations. 
The variational principle generates only a complete system of dynamic equations 
(but not their closed subsystem). When the Lin constraints 0] were added to 
the Euler equations, the system of dynamic equations became complete, and the 
variational principle was written for different partial cases ]5], || [7], ^ [| [TO], [TTL [12], 



13, 14, pl|. The variational principle was written for an ideal fluid which can be 
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considered to be a continuous set S [Sd] of nonrelativistic particles Sd, interacting 
between themselves via the pressure, which is a function of a collective variable 
(particle density) p and temperature T. Hamiltonian function of Sd has the form 

ff(t,x,p) = i£+V(x), 
2m 

The case, when the pressure is a function of p, Vp and T, was considered in the 



paper ||16|| . Now the case of arbitrary Hamiltonian function of Sd is considered. 
Parameters of the Hamiltonian (not necessarily the pressure) are considered to be 
functions of p and its space-time derivatives. 

In the second section the Hamilton variational principle is derived. In the third 
section one shows how the Hamilton variational principle generates a description in 
terms of ^-function (wave function). In the fourth section "conservative" stochastic 
systems are described on the basis of the Hamilton variational principle 



II Variational principle 

The action for the deterministic dynamic system Sd has the form 

S d : A L [x]= [ L(t,iL,dx/dt)dt, (2.1) 

A set of many independent systems Sd forms a deterministic dynamic system (sta- 
tistical ensemble) S p [Sd], whose action is a sum of actions ( |2.1| ). Let us label each 
system Sd by parameters £ = {£1, £2, •••£n}> where n is the number of generalized 
coordinates. Usually the variables £ are referred to as Lagrangian coordinates or 
generalized Lagrangian coordinates. One obtains for the action of the dynamic 
system £ p [Sd] 

„ a=n 

£ v [S d ] : A L [x] = I Po (£) L(t, x, dx/dt)dtd£, d£ = J] d£ Q , (2.2) 

a=l 

where po (£) is some nonnegative weight function. For simplicity we set po = 1, 
because in reality a use of the weight function is unessential. One can see this for a 
special case of the continuous dynamic system S [Sd] , considered in paper |TB| , where 
the weight function was used. For further calculations it is useful one to introduce 
a fictitious Lagrangian time coordinate £ — £0 (t, £) , and rewrite the action (|2.2j ) in 
the form 

/• k=n 
L(x, -)x°d n+1 Z, d" +1 e = II to, (2.3) 

x k=0 

where x = {x°,x} = {t,x} = {x k }, k = 0, l,...n, ( = {£ , = {6}, k = 
0, 1, ...n, and the point means differentiation with respect to £0 

x=^, x ^^ (2.4) 
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Now the action (|2.3|) is considered to be a functional of n dependent variables 
x. The variable t = x° is some fixed function of £. The form of this function is 
unessential. 

Let us consider the Lagrangian coordinates £ to be dependent variables, which are 
functions of independent variables - coordinates x = {£,x} = {a; ,x} = {x 1 } , i = 
0,1,. ..n and transform the action to the form, where Eulerian coordinates x are 
independent variables. Let us introduce designations 

J = d(^]t\'.'.% = det lMl> = = fc & i,k = 0,l,...n (2.5) 

l= dJ_ = d(x\^^...^ n ) p=1 o = ^L (26) 

and relations 

d£ d£ = Jd^x = J Y[ dx\ %r = Sr^TY = ^ a = l ^ - n - 

(2.7) 

The action ( |2.3| ) transforms to the form 

£ P [Si] : Ae[£) = J C (x,j) d n+1 x, (2.8) 



= L(x°,x, — )j 



i / 

0' 



where j = j} is a function of defined by the relations (p.6|). The action ( |2.8|) 
is considered to be a functional of n + 1 dependent variables £. In fact the variable 
£o is fictitious, and variation with respect to £o leads to an identity. 

Let us introduce new variables j = {j°, j} = , k — 0, 1, ...n, defined by the 
relation (|2.6| ) by means of the Lagrangian multipliers p = {pk} , k — 0, 1, ...n. 



£ p [S d ]: ^bK,p,j] = /^(^l-pilf-l-j d^, (2.9) 

Here and further a summation over repeated Latin super- and subindexes is made 
from to n. 

There are two ways for derivation of dynamic equations. The first way: the 
action (|2.9| ) is varied with respect to variables j,p,£, and after elimination of La- 
grangian coordinates £ the Eulerian equations fll.vj), ( |1.4Tj appear. We are interested 
in another way, when dynamic equations SAe/^^% = are integrated in the form 
( |L13j ) and after elimination of variables p, j, one obtains dynamic equations of the 
Hamilton- Jacobi type for hydrodynamic potentials £. 

Variation of the action (|2.9p with respect to leads to dynamic equations 



/ d 2 J \ 

S£ k : -dilPi ^z M =0 (2.10) 
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Being linear with respect to p, equations ( p.lOj ) can be solved in the form 

Pk = b (d k <p + g a (£)) d k U k = 0, 1, ...n (2.11) 

where bo is an arbitrary constant, g a (£) , a = 1,2, ...n are arbitrary integration 
functions of the variables After integration the fictitious variable £o stops to be 
fictitious, and <p is a new dependent variables which appears instead of £o- Using 
Jacobian technique [p76| , one can verify by means of a direct substitution of ( 2.11 



in (|2.10D that ( 2.1 If) is a solution of ( 2. lOf) . In particular, one needs the following 
identities 



d 2 J 



<9J dJ dJ dJ 



i,s,l = 0, 1, ...n. 



dJ 



d£ s ,k 
dJ 



SfJ, 



< dJ - x't 

£fc,2"H7 — = l J, 



i, k,l,s = 0, 1, ...n, 



I, s — 0, 1, ...n, 



Eliminating p k from relations ( |2.11|) and 



d k d£ s = 0, 



0,l,...n. 



one obtains the action 
Let us introduce new variables 



(2.12) 
(2.13) 
(2.14) 
(2.15) 

(2.16) 



v=j/j°, 



P = J 



Then the action Q2.16Q transforms to the form 

A E \p, v, <p, €] = / p{L(x, v) - & [(9o¥? + g a (£)d £a 
-bov^[d^ + g Q U)d^ a }}d n+1 x, 



(2.17) 





Let us now eliminate the variables v. Varying the action ( 2.17f) with respect to v 
one obtains 

5A E 



' dL(x, v) 



-b [dw + g a {t)d^ a }) =0, 



P = 1,2, ...n 



Let us use designations (|2.11|) . The dynamic equations (|2.18|) take the form 

/ dL(x, v) 



P0 



P = 1,2, ...71 



(2.18) 



(2.19) 



If the dynamic system <Sd is Hamiltonian, n equations (|2.19|) , can be resolved 
with respect to v in the form ir = dH/dpp, (3 — 1, 2, ...n, where 



if = #(*,x,p) 



L(x, v) 



(2.20) 
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is the Hamilton function for the dynamic system Sd- 
Substituting v in the action ( |2.17| ), one obtains 



A E [p, <p,t]=J p{-H (t, x, p) - b [d f + g a (t)d Q}d n+1 x, (2.21) 

where p is determined by the relation ( [2.11| ). 
Dynamic equations have the form 

5p: H(x,p) + b {do¥ + g a (ti)doZ a ] = (2.22) 

6<p: b °( d °P + d p(p]^)) =0 ( 2 ' 23 ) 

^ • n p (&° + w^' 1 ) = ' — (2 - 24) 



where p is determined by the relation (p. 11 
Let us consider a special case, when 

Then 

J°(«) = ^R a=l,2,..,», 

equations (|2.24 ) are satisfied identically, and p is a gradient 

Pfi = bodp (<p + <f> (£)) , /3 = 1, 2, ...n. (2.26) 

Equation ( |2.22| ) turns to the Hamilton- Jacobi equation for the variable $ = bo (if + 0) 

<%$ + # (x, p) = 0, p = V$ (2.27) 

In this case the dynamic system £ v [Sd] , considered to be a generalized fluid, has an 
irrotational vector field of momentum p (irrotational flow). 

The case, when condition ( |2.25| ) is not satisfied, the momentum vector field p 
is rotational (rotational flow). In this case equations (|2.24 ) may be considered to 
be some kind of generalization of the Hamilton- Jacobi equation. For a real particle 
whose Hamiltonian has the form 

if(x,p) = ^ + V(x) (2.28) 

the velocity v = p/m, and vector fields of p and v are simultaneously both rotational 
or irrotational. For a more general form of Hamiltonian the vector field of p may 
be irrotational, whereas the vector field of v is rotational, or vice versa. 
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Note that the action ( p.21| ) contains information on initial values of momenta p 
of dynamic system £ p [«Sd]. For instance, let us use standard initial conditions for 
variables (p, £ 

€(0,x)=x, p(0,x) = 0, (2.29) 

i.e. elements Sd of £ p [Sd] are labeled by their coordinates x at the moment t — 0. 
Then the integration functions g a (£) are determined from the relation 

p (0, x) = b g a (£) V£ a = b g (€) , g («) = {<f (£)} , a = 1, 2, ...n (2.30) 

Thus, one can consider that the action (|2.21|) contains complete information on 
motion of element Sd of £ p [Sd] ■ Vice versa, if one knows how any element Sd 
moves, one knows all about motion of the ensemble elements Sd except for density 
of elements Sd of £ p [Sd] , described by the quantity p. The variable p is a collective 
variable, containing information on density of elements Sd of the statistical ensemble 
£ p [Sd] ■ Information on initial value of p is to be given additionally. Dynamic 
equations (|2.22| ), ( 2.24 ), describing motion of single elements Sd of the statistical 



ensemble £ p [Sd] , do not depend on the quantity p. In general, evolution of the 
statistical ensemble is insensitive to the number N of its elements. This property of 
the statistical ensemble is described by the relation 

Ae[o>P, v, (p, £] = ciAe[p, v, ip, £], a = const > (2.31) 

The property (|2.31| ) is valid for any statistical ensemble £ p [Sd] and £ p [S st ]. 



Ill Description in terms of ^-function 

Let us introduce fc-component complex function ip = {ip a }, a = 1,2, ... k, defining 
it by the relations 

i> a = yfpe itp u a (i), C = v^Xte), a = 1, 2, ... k 

a=l 

where (*) means the complex conjugate, u a (£), a = 1,2, ... k are functions of only 
variables £. They satisfy the relations 

7 k Flii Pin* k 
-^E«^-^) = ^)> P = l,2,...n, 5>X = 1 (3.1) 

2 a= i C4/3 OKp a=1 

k is such a natural number that equations (|3.1|) admit a solution. In general k may 
depend on the form of the arbitrary integration functions g = {g /3 (£)}, /3 = 1, 2, ...n. 
It is easy to verify that 

p = ^*V>, ^(^a) = -^-(^^-^-^), / = 0,l,...n (3.2) 
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The variational problem with the action (|2.21|) appears to be equivalent to the 
variational problem with the action functional 



ib 



2ip*ip 



II I ./■. -— — (>*V0 - Vip* ■ ip)J ip*ip}d n+1 x 

(3.3) 

where V = {d a } , a = 1, 2, ...n 

Note that the function ip considered to be a function of independent variables 
{£, x} is very indefinite in the sense that the same state {p (t, x) , P (t, x)} of the 
statistical ensemble £ p [<Sd] may be described by different ^-functions. There are 
two reasons for such an indefmiteness. First, the functions u a (£) are not determined 
uniquely by differential equations fl3.1|). Second, their arguments | as functions of x 
are determined only to within the relabeling transformation 



Z« -> ia = D = det || d£ a /d£f, ||= 1, a, 13 = 1, 2, 



.71 



(3.4) 



Description of the statistical ensemble S p [<S<j] in terms of the function if) is more 
indefinite, than the description in terms of the hydrodynamic potentials £. Infor- 
mation about initial and boundary conditions, contained in functions g(£), is lost 
at the description in terms of the ^-function. 

Dynamic equations have the form 



sr 




ibadn - H 



dH 

-ib d - H + T—Pa 



V + V- 



dp a 



dp 




V + V- 



dp, 



dp*, 



ipp = 0, p = l,2,...k 

(3.5) 
/3 = 1,2,. ..k 







(3.6) 

where H = H (x, p) and are considered to be operators of multiplication by these 
quantities, and one has to substitute the expression ( |3.2j ) instead of p before action 
of operator V. In general, dynamic equations ( |3.5|) (|3.6j ) are not linear with respect 
to ^-function, although they may be linear in some cases. In these interesting cases 
the dynamic equations may be solved rather simply. 

The number k of the ^-function components in the actions (|3.3| ) is arbitrary. A 
formal variation of the action with respect to tp a and r/>*, a = 1,2, ... k leads to 
2k real dynamic equations, but not all of them are independent. There are such 
combinations of variations Sip a , 5ipa, a = 1,2, ... k which do not change expressions 
( 3-2|) . Such combinations of variations 5ip a , #0*, a = 1,2, ...k do not change 
the action ([3.3D, and corresponding combinations of dynamic equations 8A/8ip a = 
0, 5A/5ip* a = are identities that associates with a correlation between dynamic 
equations. Thus, increasing the number k, one increases the number of dynamic 
equations, but the number of independent dynamic equations remains the same. 
The number k is restricted underside by the condition that equations ( |3.1| ) have a 
solution. In other words, the minimal number k m of the -0-function components 
depends on the form of functions g(|), i.e. on initial conditions. This number k m 
associates with the kinematic spin ( fc-spin) s = 2k m + 1 of the ensemble state [16 
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The ^-function and fc-spin remind quantum mechanical wave function and spin 
of a particle respectively. The ^-function coincides with the wave function, if the 
dynamic equations (|3.5| ), (fOP become linear. It appears to be possible for statis- 
tical ensemble £ p [S st ] of stochastic systems «S s t. In this case the fc-spin associates 
with the spin of the described particle, but the fc-spin remains to be a property of 
the statistical ensemble £ p [S st ] ( i.e. a collective property), whereas in quantum 
mechanics the spin is considered to be a property of an individual particle. 

IV Pure Statistical Ensemble of stochastic parti- 
cles 

Let us consider statistical ensemble £ p [S st ] of stochastic particles S st . There are 
no dynamic equations for «S st , and one cannot derive dynamic equations for £ p [Sst] 
from dynamic equations for S st . But we believe that dynamic equations for £ p [«S st ] 
exist, because experiments with statistical ensembles of stochastic particles <S st are 
reproducible. Let us consider motion of «S st as a result of interaction between some 
deterministic particle S& and some stochastic agent, which makes motion of to be 
stochastic. To derive dynamic equations for E p [<S st ] , some supposition on properties 
of this stochastic agent are to be made, because one cannot derive dynamic equations 
for £ p [S s t] from nothing. If S st is a Brownian particle, moving in a gas, one supposes 
that the Brownian particle collides with gas molecules, and these collisions make the 
Brownian particle motion to be stochastic. One supposes that these collisions are 
accidental and independent. As a result a motion of the Brownian particle may 
be considered to be a Markovian process. Dynamic system £ p [<S s t] appears to be 
dissipative in this case, and there is no variational principle for it. 

However, there is another kind of stochastic agent which remains the statisti- 
cal ensemble £ p [<S st ] to be a conservative deterministic dynamic system. Dynamic 
equations of this system can be derived from the Hamilton variational principle for- 
mulated in a proper way. This stochastic agent is known as a quantum stochasticity, 
which is an origin of quantum effects, described by quantum mechanics. Stochas- 
tic systems S s t, associated with such a kind of stochasticity will be referred to as 
conservative stochastic systems. 

Let us consider some microparticle, for instance, electron. An individual electron 
is a stochastic system «S st , because experiments with a single electron are irrepro- 
ducible. For instance, in a diffraction experiment a single electron, flying through a 
narrow hole in a diaphragm, hits a new point of a screen any time. But distribution 
of many independent electrons over the screen surface is reproducible. It means, 
that the statistical ensemble £ p [S st ] of many <S st is a deterministic dynamic system 
Ss- This deterministic dynamic system <Sg is described by the action 




x 
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where ip is so called wave function. Dynamic equation, generated by the action 



h 2 

ihd ip + — V 2 V» = (4.2) 
2m 

is known as the Schrodinger equation. 

Why an electron is described in terms of wave function? What is the wave 
function? Principles of quantum mechanics answer these questions. The wave func- 
tion is a fundamental mathematical object (a point in the Hilbert space) which 
evolves according to linear dynamic equation. All relations of quantum mechanics 
are obtained from corresponding relations of classical ones, replacing position x — ► 
operator of multiplication by x and momentum p — > — ihV. Quantum mechanics is 
an axiomatic theory, and quantum principles work very well in the non-relativistic 
case. Attempts of extension of the quantum mechanic principles to the relativis- 
tic case met problems. Investigators were forced to introduce corrections and new 
principles. All this leads to a suspicion that quantum mechanical principles are con- 
ceptually non-relativistic, and principles of quantum mechanics cannot be combined 
with the relativity principles. 

We shall not use quantum mechanical principles. Instead we suppose, that the 
quantum stochasticity has a geometrical origin. Namely, we suppose, that the real 
space-time is described by the geometry of Minkowski only approximately. This ap- 
proximation is valid only for large enough values of world function for the Minkowski 
space-time 

a M {x, x) = a M (t, x, t', x) = - (c 2 (t - t'f - (x - x ) 2 ) (4.3) 

where x = {t, x} and x' = {t', x'} are coordinates of two points in the space-time, 
and c is the speed of light. The world function a is a way of geometry description 
17] . Recently it was shown that the world function describes geometry completely 



fig , fL9|l , and any change of world function leads to a change of geometry and vice 
versa. Geometry of the real space-time is described by the world function 

cr (x, x ) = a M (x, x') + D (a M (x, x')) (4.4) 
Here <tm is the Minkowski world function (|4.3j), D is a distortion function 



D= \d if a M > op = h = ^ l - 17 g/cm (4.5) 

\ if a M < 2bc &/ v ; 

where H is the quantum constant, and b is a new universal constant. <jq is a constant 
of the order d ~ 10~ 21 cm 2 . Correction to the Minkowski world function is small 
for large values of aM, because the characteristic length \fd ~ 10~ n cm is essential 
only in the microcosm. Geometry, described by the world function ( f4.4|) , is non- 
Riemannian. (It is called tubular geometry, or briefly T-geometry). In T-geometry 
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a motion of free particles is stochastic, although the world function (fO) and T 



geometry in itself is not stochastic^]. The stochasticity in T-geometry, described by 
the world function Q4.4Q , depends on the particle mass. The stochasticity is large 
for particles of small masses. It is negligible for particles of macroscopic mass. 

Construction of a geometry on the basis of only world function is a new concep- 
tion of geometry. (See for details Ref. |1| ^(J, ^J, |!9lD. Now from this conception 
we need only existence of a geometrical stochasticity. We identify the geometrical 
stochasticity with the quantum stochasticity and try to construct a statistical de- 
scription of microparticles on the basis of Hamiltonian variational principle. Such 
an approach has many points in its favour. First, neither quantum principles, nor 
new principles are used. All results are logical corollaries of supposition about the 
form of the world function ( |4.4|) of the space-time. Second, such an approach is rela- 
tivistic from the outset, and one does not need to combine relativity principles with 
quantum mechanical ones and to solve problems, connected with this integration. 
Third, the wave function and spin appear to be attributes of statistical descrip- 
tion (not fundamental objects, whose meaning is obscure). Fourth, interpretation 
of quantum mechanics appears to be single-valued. Different versions of quantum 
mechanics interpretation, which take place in the axiomatic quantum mechanics, 
are impossible under such an approach. In particular, there is a clear distinction 
between the individual stochastic system S st and statistically averaged system (<S st ) . 
(The last is the statistical ensemble £ p [S st ] normalized to one system). In quantum 
mechanics the same term is used for both S st and (S st ). It leads to misunderstand- 
ing and paradoxes such as paradox of the Schrodinger cat, or EPR paradox. Fifth, 
a use of the model approach (instead of axiomatic one) admits one to determine 
the domain, where the quantum principles are valid, and quantum mechanics works 
correctly. 

Let now S s t be an electron or some other microparticle. The system <S s t is stochas- 
tic, but statistical ensemble £ p [S st ] of electrons S st is a conservative deterministic 
dynamic system, and some variational principle takes place for £ p [S st \ . We suppose 
that in the zero approximation £ p [<S st ] is described as a statistical ensemble £ p [S^] , 
where is a deterministic relativistic particle, described by the Hamiltonian 



H (x, p) = ^m 2 c A + p 2 c 2 (4.6) 
Then variational principle ( p.21| ) for dynamic system £ v [<Sd] has the form 



A E \p,<P,t] = J p{-y/mW + p 2 c 2 - b [d ip + g a U)d £ a }}d 4 x, (4.7) 



where p is determined by (|2.11| ), n = 3. To obtain variational principle for £ p [S st ], 



one should take into account stochastic component of the particle momentum. The 
mean value p s t of the stochastic component has the form 

p st = Wlnp (4.8) 



3 It seems rather evident that the free particle motion is stochstic in the space-time with stochas- 
tic geometry |2^, |2^, |24| , but a stochastic motion of a free particle in the deterministic space-time 
looks rather unexpected. 
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where % is the Planck quantum constant. In other words, the mean value p st of the 
stochastic component depends on the collective variable p which describes the state 
of the statistical ensemble. As far as the stochastic components and regular ones are 
supposed to be independent, a sum of squares of these components is to be used in 
the modified Hamiltonian. 

The action for the dynamic system S p [S st ] has the form 

A E \p,<p,t] = J p{-\lmW + p 2 c 2 + h 2 c 2 (V\npf-b [d oV + g a (i)do^}}d 4 x, (4.9) 



where p is determined by ([2.11 ) , n = 3. This fact may be interpreted also in the 



sense that the mass of the electron is modified 

n 



2 



m 2 ^ m 2= m 2 + — (Vlnp) 2 (4.10) 

q C 



In fact, the supposition (|4.10|) is an origin of supposition (|4.8|) . One supposes that 



consideration of geometric stochasticity leads only to a change of parameters of the 
system iS d in £ p [«S d ] . These parameters start to depend on the state density p of the 
ensemble. The system <S<j ( |4.6|) has only one parameter - a mass m. The dependence 
( |4.10| ) is obtained as a result of averaging over world lines of stochastic particles p0 



It is to be invariant with respect to a change of number of the ensemble elements, 
i.e. with respect to transformation 

p — > ap, a = const > 0. (4-H) 

Note that the action ( |4.9| ) cannot be considered to be an action for the statistical 
ensemble of any deterministic dynamic systems <S<i, because the effective Hamiltonian 



H eS = ymV + p 2 c 2 + h 2 c 2 (Vlnp) 2 , (4.12) 

which enters in the Hamilton variational principle, depends now on the collective 
variable p describing the state of the whole ensemble, whereas Hamiltonian for any 
statistical ensemble £ v of deterministic dynamic systems has to depend only 
on variables of the dynamic system as it follows from ( 2.21| ). The action (|4.9|) is 



an action for a set of identical deterministic dynamic systems Sd, interacting between 
themselves. Hence, it is not a statistical ensemble of S^- In the same time the action 
( O ) may be considered to be an action for a statistical ensemble, because it has the 



main property of a statistical ensemble: not to depend on the number of elements 
of the statistical ensemble. Mathematically it means that if p is substituted by ap, 
a =const, Ae[p, <p, £] is substituted by <iAe[p, <p, £]> or 

A E [ap, <p, £] = clAe [p, V, €] (4-13) 



It is easy to see that the action ( [4.9|) has the property (|4 .13 ). Hence, the action 



( |4.9| ) is an action of a statistical ensemble, but this ensemble cannot be a statistical 
ensemble of deterministic dynamic systems, because its elements interact between 
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themselves and are not independent. It means that the statistical ensemble ( f4.9|) is a 
statistical ensemble of nondeterministic dynamic systems. In other words, a statisti- 
cal ensemble of nondeterministic systems can be considered to be a set of interacting 
deterministic systems, i.e. a stochasticity of dynamic systems is substituted by in- 
teraction of deterministic dynamic systems. Form of this interaction depends on 
the form of stochasticity. Considering different forms of interaction, satisfying the 
constraint ( 4.13Q , one can label and investigate different forms of stochasticity. 



In fact we have no other way of the stochasticity description except for a reduc- 
tion of the statistical ensemble £ p [S st ] to a set iS re d [Sd] of interacting deterministic 
dynamic systems Sd- 

£ p [S st ]=S rcd [S d ] (4.14) 

Character of stochasticity is described by the way of interaction between determinis- 
tic dynamic systems Sd in the set S Ted [Sd] which is a deterministic dynamic system, 
consisting of interacting deterministic systems Sd- We shall refer to S ved [S d ] as a 
reduced statistical ensemble, consisting of elements S d , although the set S red [S d ] of 
S d is not a statistical ensemble of S d at all. 

Thus, we know the only way of effective mathematical description and investi- 
gation of stochastic dynamic systems. This is a substitution of stochasticity by an 
interaction, i.e. reduction of the statistical ensemble £ p [S st ] to a set S red [S d ] . One 
can describe properties of stochastic system S st , only referring to the properties of 
the reduced statistical ensemble iS re d [Sd]- Different kinds of stochasticity are de- 
scribed by a consideration of different types of interaction between S d in S Ted [S d ]. 
We shall show that the action ( |4.9|) describes a dynamic system which in the non- 



relativistic approximation is described by the Schrodinger equation. But there are 
another reduced statistical ensembles S rcd [S d ] which have the Schrodinger equation 
as a dynamic equation. 

For instance, the reduced statistical ensemble S ved [S d ] , described by the action 

A E \p,<P,iA = J pi-y/mW + h 2 c* + k'k { ) - b [d p + g a ($)d Z a ]}d n+1 x, 

(4.15) 

k = {k , Ki, k 2 , k 3 } , k 1 = g ^Kj, 1 = 0,1,2,3 (4-16) 
where S d interact via some relativistic quantum K-field, is also described (under 



some conditions) by the Klein-Gordon equation |25|. In the non-relativistic approxi- 
mation this equation reduces to the Schrodinger equation. There are another kinds 
of a reduced statistical ensemble <S re d [Sd] which has the Schrodinger equation as a 
dynamic equation under some conditions. It is not clear a priori, which of these 
reduced statistical ensembles S red [S d ] is true. This question needs further investi- 
gation. 

Let us return to the action Q4.9|) and represent it in terms of ^-function. In gen- 
eral, introduction of interaction between Sd in the action (|4.9|) is not relativistically 



covariant. In the non-relativistic approximation the action (|4.9|) has the form 
A E \p,<P,i\ = J P{~™ 2 - |1- |^(Vlnp) 2 -bo[d p + g a (H)d Z a ]}d 4 x, (4.17) 
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where p is determined by the relation (p. 11 



In terms of the ^-function the action (4.17) is written in the form 



r ib 0f i*a i * /* n 2 fr 2 (Vp)' 



+ 



b 2 



ipm 



(4.18) 



where p = ■0* , 0. 

Let the function t/> have components. Regrouping components of the function 
if) in the action ( 4.18] ), one obtains the action in the form 

+ T E <E* >7 <W + 4^ (Vp) 2 - mc 2 p}d 4 x, p = W (4.19) 



Q,/3=l 



6 pm 



where 



Qa, 



13,1 



4>a Ipfi 



a,P = \,2,...k 7 = 1,2,3 



(4.20) 



and Q* a p a is complex conjugate to Q a/3,7 ■ 

In the simplest case, when the ^-function has only one component, all quantities 
Qn,7 = 0, 7 = 1,2,3, and the motion of the ensemble particles is irrotational. 
Then the action (|4.19| ) reduces to the form 



mc 2 p 



\^ (v P ) 2 }d 4 x, p = w 

8pm 



(4.21) 



Dynamic equation, generated by the action ( |4.21|) is nonlinear due to the last 
term in the action ( |4.21| ). Equating the arbitrary integration constant b to h 
(bo = h) , one obtains instead of (|4.21|) 



At 



[fa fa] = I { — (ip*d ip - d ip* ■ if)) - —Vfa ■ W - mc 2 fa ■ ^}d 4 x (4.22) 



Note that the equivalence of the actions (|4.21 ) and (|4.22|) can be verified directly, 
making in ([4.21]) a change of variables 



V^^H^I exp \J- log 



(4.23) 



which transforms the action Q4.21 ) into the action Q4.22j ) 
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The action (|4.22|) generates linear dynamic equation for the ip- function. 



6 A _ „ ft 



2 



iW Q ip + — V 2 ^ - mc 2 tp = (4.24) 



Sipa 2m 



After substitution ^ = ^exp [^^J the dynamic equation (4.24) transforms to the 
Schrodinger equation for a free particle 

ft 2 

ihdoV + — V 2 ^ = (4.25) 
2m 

It means that in some cases the ^-function can coincide with the wave function. 
Let us compare actions ( 4.21| ) and ( |4.22[ ). They differ in a choice of the integra- 



tion constant b and describe the same dynamic system. The action (|4.21|) contains 
only one quantum term (i.e. the term, containing the quantum constant ft). This 
term — ft 2 (Vp) 2 /8pm describes the energy density of stochastic component of mo- 
tion. Other terms, containing integration constant bo, are usual dynamical terms. In 
the action (|4.22|) practically all terms contain the quantum constant ft, and should 
be interpreted as quantum terms. Extension of quantum properties to all dynamical 
terms is the price which is paid for linearity of dynamic equation ( |4.24j ). 



Separation of quantum and classical properties is very simple in the action ( |4.21|) 



Setting ft = in (|4.21), one suppresses all quantum properties, and remaining terms 



describe a pure ensemble of deterministic classical particles. In the action (|4.22 ) one 



cannot set ft = 0, because all dynamic terms contain ft. Separation of the classi- 
cal part of the action in ( f4.22|) is rather complicated procedure. From viewpoint of 
statistical description a linearity of dynamic equation (|4.24j) is an accidental circum- 
stance, connected with a special form — ft 2 (Vp) 2 /8pm of the stochastic component 
energy density and with irrotational character of the quantum fluid flow. If this flow 
is rotational, and ^-function has more, than one component, dynamic equations are 
not linear. One cannot be sure, that completely relativistic statistical description 
can reduce to a linear dynamic equation, because in this case the stochastic com- 
ponent energy density has the form, which differs from — ft 2 (Vp) 2 /8pm. Neverthe- 
less there exists such a relativistic statistical description, which generates a linear 
dynamic equation for one-component ^-function (the action ( |4.15| ) generates the 
Klein-Gordon equation for one-component ^-function ||25|| ). 

Note that the term pU st = ft (V p) 2 /8pm, describing stochastic component energy 
density, becomes to be large, if there is a local increase of density p. Let 5x be 
the linear size of the region of the density increase, then |Vp| ~ p/Sx, and the 

H 2 

stochastic component energy U st ~ &m ^ x f becomes very large, provided 5x is small 
enough. After expansion of the quantum fluid the local density increase disappears 
and the stochastic component energy U st transforms to the kinetic energy p 2 /2m. 
The particle momentum p becomes to be of the order p « V 2mU st ~ This 
relation is an origin of the indeterminacy relation. It means that the well known 
indeterminacy relations is a corollary of the way as the stochastic component energy 
U st depends on the statistical ensemble density p. 
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Of course, linearity of dynamic equations is very convenient and useful, but 
it does not mean that linearity of dynamic equations may be considered to be a 
principle for construction of relativistic quantum theory. 
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